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During the 1996 HST/FOS monitoring campaign of the Seyfert 1 galaxy NGC 3516, 
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the UV continuum showed a factor of 5 variation; correspondingly the strongest broad UV 

h ; 

emission lines (e.g. Lya, C iv) displayed somewhat weaker, though still significant variations 
(a factor of 2). In contrast, no variation was detected in the strong broad Mgll + UV Fell 

•rH . 

^ ■ emissiondine complex. While historically the flux in the broad Mg II emission line has varied 
(a factor of 2 within a few years), the profile shape has not been observed to change over the 
past decade. In contrast, the high ionization lines (HILs) show both emission-line flux and 
profile changes on relatively short timescales, which appear to be correlated with changes 
in the continuum level. Significantly, on 1996, February 21, the profile shapes of the HILs 
(e.g. Lya, Civ) and low ionization lines (LILs, e.g. Mgll) were identical, at which time the 
continuum level was at its highest. These results impose stringent constraints upon viable 
kinematic models of the BLR in NGC 3516, some of which we discuss below. 

Subject headings : galaxies : Seyfert — galaxies : individual (NGC 3516) — galaxies : 
emission lines — galaxies : kinematics and dynamics 
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1. Introduction 

NGC 3516 is a Seyfert 1 galaxy with a well documented history of short-timescale, 
large-amplitude continuum and emission-line variability (Koratkar et al. 1996). Historically, 
IUE/SWP spectra of NGC 3516 have also displayed what has previously been referred to as 
a strong, variable, broad absorption line (VAL, FWHM ~200 km s" 1 ) (Kolman et al. 1993; 
Walter et al. 1990; Voit et al. 1987), presumed to be similar to that observed in Broad 
Absorption Line quasars (BALs, Weymann et al. 1991). The subsequent disappearance of 
the supposed VAL between 1989 and the onset of the 1993 IUE/SWP monitoring campaign 
(Koratkar et al. 1996), may be associated with the increased ionization state of the gas 
inferred from the observed variation in the strength of the X-ray warm absorber detected by 
ROSAT (Mathur et al. 1997). 

However, recent high resolution UV observations of the Civ emission-line region in 
NGC 3516 with the HST/GHRS (Crenshaw et al. 1998), show that NGC 3516 displays 
multiple intrinsic absorption lines that are all blue-shifted and relatively narrow (FWHM 
~500 km s^ 1 ), and which appear to remain invariant on timescales of ~ 6 months, and 
are thought to originate in the narrow-line region gas. In light of these observations, Goad 
et al. (1999) argue that the supposed VAL is in fact multiple, narrow absorption- line systems 
which appear broad at the lower resolution of IUE. In this picture, the historical variation in 
the strength of the VAL results from variations in the number of narrow-line clouds coupled 
with changes in the shape of the underlying broad emission-line. 

In an ongoing effort to determine the nature of the broad emission/absorption line regions 
and their long-term evolution, we obtained 5 HST/FOS cycle 6 observations of NGC 3516, 
from 1995 December to 1996 November, covering the wavelength range AA1 150-3300 A. In 
this paper we present the first results of this study, highlighting the absence of emission-line 
flux variations in the broad Mg II and UV Fe II emission lines, and its implications in terms 
of viable kinematic models of the BLR. 
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1.1. Results from the 1996 HST/FOS monitoring campaign 

Goad et al. (1999) showed that during the 1996 HST/FOS monitoring campaign, 
NGC 3516 displayed the following continuum and emission-line variability characteristics: (i) 
large amplitude, wavelength dependent, continuum flux variations (a factor of 5 at A1365 A, 
c.f. a factor of 3 at A2230 A), (ii) correlated, large amplitude, UV broad emission- line flux 
variations (a factor of 2 for both Lya and C iv) , and (iii) the complete absence of variations in 
both the strength and shape of the broad Mg II and UV Fe II emission-line complex (Figure 1). 

A comparison of the Mgll emission-line strength determined from IUE/LWP observations 
taken in 1993 (Koratkar et al. 1996) with the present data, show that the flux in Mgll has 
apparently remained constant for at least 4 years. The absence of significant variability in the 
Mg II emission-line strength, despite large amplitude continuum and HIL flux variations, has 
also been noted in lower resolution (pre-1989) IUE/LWP data (Koratkar et al. 1996) when 
the Mg II emission-line strength was a factor of 2 weaker. Interestingly, the observed doubling 
of the strength of the broad Mgll emission-line between 1989-1993, appears to coincide with 
the disappearance of the higher-velocity narrow absorption-line components comprising the 
so-called VAL. It remains to be seen whether these two phenomena are related. Significantly, 
despite the notable increase in the Mgll emission-line strength between 1989 and 1993 its 
shape has remained unchanged for over a decade (Figure 2). Moreover, in the present 
campaign, the shapes of the broad Lya, Civ and Mgll emission lines as observed on 1996, 
February 21, when the continuum was in its highest state, are indistinguishable from one 
another (Figure 3a, N.B. Lya is not shown). 



2. Implications for kinematic models of the BLR 

A viable kinematic model of the BLR in NGC 3516, must explain : (i) the absence 
of both short-timescale (~ months) and long-timescale (~ a few years) variations in the 
Mg II emission-line flux and shape, despite significant changes in the strength of the ionizing 
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continuum; and (ii) the remarkable similarity between the line shapes of the HILs and 
LILs, at high continuum levels, when clearly these lines form under very different physical 
conditions. 

2.1. Flux variability constraints 

An absence of response in the Mgll emission line can result for a variety of physical 
reasons : (i) the continuum band driving the line is invariant; (ii) the emission-line is 
insensitive to continuum variations; and (iii) the line emitting region is physically extended. 

Here we address the shortcomings of each of these scenarios. 

(i) While it is true that Mgll and Lya respond to very different continuum bands (Mgll is 
mostly sensitive to the continuum from 600-800 eV (Krolik and Kallman 1988) whereas Lya 
is driven mainly by Lyman continuum photons), simultaneous multi- wavelength observations 
of a small number of AGN in the UV, optical and soft X-ray band-passes, indicate that 
the amplitude of the continuum variations generally increase toward shorter wavelengths 
(Romano and Peterson 1998). Although we cannot confirm whether such a trend holds at 
600-800 eV, the Seyfert 1 galaxies for which simultaneous observations were made in the 
UV/EUV, NGC 5548 (Marshall et al. 1997) and NGC 4051 (Cagnoni et al. 1998), showed 
increased variability at shorter wavelengths. Thus, the evidence suggests that the lack of 
Mg II and Fe II variations is not due to an absence of continuum variability in the bandpass 
responsible for driving these lines. 

(ii) Goad et al. (1993, 1995) and O'Brien et al. (1995) demonstrated that the line 
responsivity r], the fractional change in line emissivity for a given fractional change in the 
continuum level, is relatively modest for Mgll (t] ~ 0.2 — 0.3), when compared to the 
approximately linear response (7/ ~ 1) of the HILs. This differences arises because for a 
single cloud and a fixed continuum shape, the variation in the position of the hydrogen 
ionization front is generally small compared to the overall extent of the partially ionized 
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zone where LILs such as Mg II are formed. For a fixed continuum shape this condition holds 
true over a large range in ionizing continuum luminosity (factors of 10 or more). Although 
the Mg II line responsivity is small, given the factor of 2 change in the C IV emission-line flux, 
we should have detected ~ 40% variation in the Mgll emission-line flux. Since no significant 
(<7%) variation was detected, low emission-line responsivity cannot be the sole explanation 
for the absence of Mg II emission-line variations. 

(iii) Finally, the lack of response in the Mg II emission line can also result if the Mg II line- 
emitting region is physically large, and has thus yet to respond to the observed continuum 
variations. While plausible, given the short-timescale response of the C IV emission-line to 
the continuum variations (~ 4.5 days, Koratkar et al. 1996), it is difficult in this picture to 
explain the similarity in shape of the C IV and Mg II emission-line profiles at high continuum 
levels, particularly if the size of the regions where these lines originate differ by a few orders 
of magnitude. 

The similar ionization potential of Fe II and Mg II suggests that these lines are formed in 
similar regions, hence it is unsurprising that the Fe II emission lines display a similar lack of 
response to continuum variations (N.B. relatively few of the Fell lines (although strong) are 
resonance lines, although others are pumped by the continuum, especially in the presence 
of extra thermal particle speeds). The smoothness of the Fell complex in NGC 3516 when 
compared to similar observations of I Zw I (Goad et al. 1999), a quasar with strong narrow 
Fell emission-lines (Laor et al. 1997), supports this finding, and suggests that the Fell lines 
in NGC 3516 are indeed produced in high velocity gas (FWHM~ 4500 km s _1 ). 



3. Profile variability constraints 

The similarity in shape between the HILs and LILs at the highest continuum levels 
suggests that in the high-state both Civ and Mgll arise in kinematically similar regions. 
However, if the lack of response of Mg II is not the result of an absence of variability in the 
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continuum band driving this line, then the difference in response timescale between the HILs 
and Mg II suggest that these lines do arise in spatially distinct regions, with Mg II formed in 
a region which extends over several light-months or more. Figure 3b indicates that between 
the low and high-states, the core of the C IV emission line shows a deficit of response redward 
of line center from to +4000 km s -1 (even after allowing for possible contamination by the 
variable absorption-line and narrow emission-line components) and an enhanced response in 
the far red-wing (> 5000 km s" 1 ). This may result from : (i) reverberation within a radial 
flow; (ii) a change in the radiation pattern of the ionizing continuum; or (iii) a change in the 
spatial distribution of the line emitting gas. 

Taken together, this evidence places severe constraints upon the spatial distribution and 
kinematics of the line-emitting gas in NGC 3516. For example, based on this data the 
following models can be excluded : 

(i) Radial flows at constant velocity, illuminated by either an isotropic or anisotropic 
continuum source (Goad and Knigge 1999). 

(ii) Radial flows in which Mg II and C IV arise in azimuthally separated cloud populations, 
either through spatial segregation or anisotropic illumination. 

(iii) A spherical distribution of clouds in Keplerian motion and illuminated by an isotropic 
continuum source. 

We propose two models which may account for the observed emission-line variations: (i) 
an anisotropic continuum source model and (ii) a terminal wind model (i.e. a flow that 
accelerates from rest until it reaches constant velocity). 

(i) The anisotropic continuum model of Wanders et al. (1995) can broadly reproduce the 
gross details of the observations reported here. In this model, the continuum is assumed to 
be comprised of two components, an anisotropic variable component responsible for driving 
the variations in the HILs and an isotropic non-variable component responsible for producing 
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the Mg II and Fe II emission. To match the observed variations, the BLR must be comprised 
of two distinct cloud populations, a low column density population producing significant 
HIL emission, and a spatially distinct high-column density population producing both HILs 
and LILs. The variable ionizing continuum component preferentially illuminates the low- 
column density population, giving rise to the variations in Lya and Civ emission-lines, 
whereas the LILs arise predominantly from the high-column density clouds illuminated by 
the non- variable isotropic continuum component. If the ionization parameter U is large 
enough, and the EUV is not absorbed, these high column density clouds could also produce 
significant HIL emission. To account for the similarity in profile shape in the high-state, 
we further assume that the variable ionizing continuum component becomes more isotropic 
when brighter, so that the C IV emission arises predominantly from the high-column density 
clouds. Since the FWHM of the low and high-state Civ emission-line profile is effectively 
the same, the velocity field cannot be predominantly radial, otherwise a different range in 
projected velocity for the HILs and LILs will result. Instead a randomized or chaotic flow 
is favored. While plausible, this model requires a high-degree of fine tuning to produce 
the detailed profile changes observed in the HILs. Hardest to explain are the reversals in 
HIL profile shape, from red-asymmetric at low continuum levels to blue-asymmetric at high 
continuum levels (Figure 3b). 

(ii) One possible scenario which may account for the above observations is a terminal wind 
model for the BLRQ. A terminal wind is the one physical structure which can display the 
same range in projected velocities on both small and large spatial scales. The wind may be 
spherical or non-spherical. For a non-spherical wind, for example a bi-conical flow, similar 
profiles for the HILs and LILs will result regardless of their respective radial emissivity 
distributions, provided that the opening angle of the cone remains constant. However, to 

x The best available variability data on NGC 3516 (Wanders et al. 1993; Koratkar et al. 
1996) suggest that the BLR kinematics are not dominated by radial motion, but may be 
consistent with a mixture of velocity components. 
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produce profiles similar to that observed here (i.e. broad wings and narrow cores), we require 
a strong azimuthal dependence on the velocity field, or a significant circularized velocity 
component to the flow (Goad and Knigge 1999). That is, the simplest terminal flow model 
in which the velocity field is constant everywhere cannot produce the correct profile shapes. 
However, hydromagnetic wind models (e.g. Emmering et al. 1992; Bottorff et al. 1997) 
do exhibit these basic properties. While the emission-line flux variations of the HILs are 
assumed to be driven by continuum variations, we propose that in this model the changes in 
profile shape at low continuum levels are in part due to changes in the formation radius of 
the HILs. Specifically, we propose that in the low-state, conditions within the BLR gas are 
such that a large fraction of the C IV emission arises at the base of the wind, in the region 
where the gas begins to accelerate (Goad and Knigge 1999). 



4. Conclusions 

The picture we envisage is one in which the steady-state profile is dominated by the 
dynamics of the system, with lines formed in kinematically similar regions, possibly a 
terminal flow. The stability of the Mgll profile on timescales of several years, provide 
strong supporting evidence for the existence of such a structure. Moreover, the absence of 
significant flux variations in the Mgll line on timescales of a few years suggest that this 
region is physically extended in size. The historic variation in the Mgll emission-line flux, 
but not its profile, could be the result of an accretion event, a disk instability, or an increase 
in the wind density and may be linked to the disappearance of the Civ VAL. 
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Fig. 1. — The high-state (epoch 2, solid line) and low-state (epoch 5, dotted line) continuum 
subtracted Mgll+Fell emission-line complex and their difference (dashed line). For clarity, 
each spectrum has been binned into 2 A wide bins. The derived errors (poissonian) in the 
difference spectrum, displaced for clarity, are indicated by the vertical bars. We detect no 
significant variation in the broad Mg II and UV Fe II emission-line complex. The absorption 
lines are Galactic in origin. 

Fig. 2. — A comparison between the normalized 1996 HST/FOS Mgll emission-line profile 
(smoothed to the resolution of the IUE/LWP), with that determined from a) the 1993 
IUE/LWP monitoring campaign, and b) the IUE/LWP archival data. Note that for the 
upper panel no scaling between the HST and IUE spectra has been applied, the emission- 
line fluxes and profile shapes are the same. The solid line in the lower panel shows the 
old average IUE/LWP Mgll flux scaled relative to the flux measured with HST. While the 
Mgll emission-line flux clearly varies (a factor of 2.3), the profile shape remains invariant 
on timescales of several years. The absorption dip blueward of line centre is due to Galactic 
Mgll. 

Fig. 3. — (a) A comparison between the Civ (solid line) and Mgll (dashed line) emission- 
line profiles as observed on 1996, February 21. Each line has been transformed into velocity 
space, and normalized to the total flux in a 2000 km s _1 wide bin centered at +2000 km s _1 . 
Note the good agreement between the line profile shapes in the blue and red wings of the 
lines. The far blue wing of Mg II is contaminated by blended Fe II hence the poor agreement 
in this region, (b) The high-state (epoch 2, solid line) and low-state (epoch 5, dotted line) 
continuum subtracted C iV+He II emission lines, scaled to the integrated flux in a 2000 km s _1 
wide bin centered at —4000 km s -1 . The most significant shape changes occur redward of 
line center, with the core displaying a deficit of response, whereas the far red wing displays 
an enhanced response. A similar result is found for Ly<x 
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